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We present a comparison of three different types of Langevin equation 

exhibiting absorbing states: the Langevin equation defining the Reggeon field 

ON 

l/~) ' theory, one with multiplicative noise, and a third type in which the noise 



is complex. Each one is found to describe a different underlying physical 
mechanism; in particular, the nature of the different absorbing states depends 






& ' on the type of noise considered. By studying the stationary single-site effective 

a" 

potential, we analyze the impossibility of finding a reaction-diffusion model 



o 
o 



ctf 



in the multiplicative noise universality class. We also discuss some theoretical 
questions related to the nature of complex noise, as for example, whether it is 



/\ • necessary or not to consider a complex equation in order to describe processes 



as the annihilation reaction, A + A — ► 0. 

PACS: 

Different systems and models appearing in physics as well as in other fields can exhibit 
absorbing states. An absorbing configuration is one in which a system can get trapped, 
from which it cannot escape |T|-§[. Therefore an absorbing configuration is a fluctuation- 
free microscopic state. Some examples of systems exhibiting absorbing states, among many 
others, are: chemical reaction-diffusion models of catalysis [|j], models for the spreading of 
epidemics or forest fires 0, directed percolation PH, the contact procces |7|,||], models of 
branching and annihilating random walks ||, damage spreading HI, and even self-organized 



systems [|IIJ . 



As some control parameter is changed, many of these systems experience a phase transi- 
tion from an absorbing phase, i.e., a phase in which the absorbing state is the only stationary 
state 0, to an active phase, characterized by a non- vanishing value of the order parameter. 
At the critical point, these systems exhibit universal features. 

It was conjectured some time ago by Janssen and Grassberger |I1J that all the different 



systems and models with an unique absorbing state, a single-component order parameter, 
and no extra symmetry or conservation law belong in the same universality class as directed 
percolation (DP), which is considered the canonical representative of that vast class of 
models. In a field theoretical description this universality class is represented by the Reggeon 
Field theory (RFT), [12) which in terms of a Langevin equation reads [p~3|j : 

a' = V2ra ( x ' *) + ara ( x ' *) ~ bn2 ( x ' t) + \f n (x, t)r](x, t) , (1) 

where n(x, t) is a density field at position x and time t, a and b are control parameters, and 
i](x,t) is a Gaussian noise which only non-vanishing correlations are: < rj(x,t)rj(x',t') >= 
D5(x — x')5(t — t') . The equation is interpreted in the Ito sense [TJ|]. The noise term in eq.(|T]) 
is proportional to the square root of the field, therefore in the absorbing state n(x) = the 
dynamics is completely frozen: both the deterministic and the stochastic terms are equal to 
zero. Other higher-order terms could be added to the deterministic part of eq. (H) but they 
can be easily argued to be irrelevant in a renormalization group sense. 

The previous conjecture has been confirmed in a large number of computer simulations 
and series expansion analysis, and DP-universality class has proven to be extremely robust 
against the modification of many details in the microscopic models. The conjecture of 
universality has been extended for multicomponent systems [15] , as well as for systems with 
an infinite number of absorbing states [|16j . 

Nevertheless, not all the systems with absorbing states belong in the universality class 
of DP. Some other classes different of DP have been identified, all of them showing some 
essential physical differences with DP. 

Two relevant examples for what follows are: 



1. Particle systems in which evolution occurs only at the interfaces separating occu- 
pied (active) from empty (absorbing) regions belong to the compact directed percola- 
tion (CDP) universality class ||17|| . Examples of this universality class are the one- 
dimensional diffusion-limited reactions of the type A + A — > and A + A — > A (in 
dimensions larger than one these models are not expected to be in the same univer- 
sality class as CDP). The pseudoparticles A can be thought of as the kinks separating 
active from inactive regions, where the dynamics occurs. A field theoretical descrip- 



tion for such class of systems was proposed by Peliti in jl8| , |19| ; its equivalent Langevin 
equation reads: 

dn(x, t) 



at 



V n(x,t) — bn (x,t) + in(x,t)rj(x,t) , (2) 



where n(x,t) is a field, % is the imaginary unit, and r\ is a Gaussian noise with some 
amplitude D. Note that n(x) = is an absorbing configuration. Hereafter we refer 
to eq. (0) as Peliti's field theory. It is important to point out that the field n(x, t) is 



not the density of A particles, but a more abstract field |20 which expectation value 
coincides with that of the real density field, and which higher order moments can be 
also related to higher order moments of the density field. All this will become clearer 
in a forthcoming section where the explicit derivation of eq. @ is performed. 

Before concluding this epigraph let us point out that there is a class of systems which 
present an active as well as an absorbing phase with a phase transition separating both 
of them, and which associated noise should also present a complex structure: this is 
the so called parity-conserving universality class ||. Due to the extra conservation law 
it is clear by now that these systems are not in the RFT universality class. 

2. A new universality class characterized by a noise different from that of the previously 
described classes has recently been elucidated ( see |j2T^f23| and references therein) : the 



multiplicative noise (MN) universality class. While in the Langevin equation for the 
RFT the noise amplitude is proportional to the square root of the field at each point, 
in the MN universality class the noise amplitude is proportional to the field itself, i.e. 



dn{x,t) _„ L ,, ,., ...,,,.,., , , 



ctt 



V n(x, t) + an(x, t) — bn (x, t) + n(x, t)r)(x, t) , (3) 



with n(x, t) being a density field at position x and time t, and r](x, t) a Gaussian white 
noise. Obviously this new type of noise is also compatible with the presence of an 
absorbing state at n(x) = at which the dynamics is completely suppressed. 

While there are many microscopic react ion- diffusion models belonging in the DP univer- 
sality class, and is also easy to identify reaction-diffusion models in the Peliti's field theory, 
no microscopic reaction-diffusion model in the MN universality class has been identified so 
far. The possibility of constructing a reaction-diffusion model that exhibits the rather strik- 



ing properties of MN [22] has been explored in a recent paper by Howard and Tauber [24 
They concluded that given the apparent impossibility of finding such type of model in the 
MN class, the physical meaning of that universality class is unclear. 



Motivated by the previous work [24| we have further investigate this issue. In what follows 
we present a comparison of the different noise terms appearing in Langevin equations for 
Reggeon field theory, the multiplicative noise, and Peliti's field theory, to identify physical 
differences among them. A simple and intuitive justification of the fact that no reaction- 
diffusion system can be found in the MN universality class is given. Alternatively, we 
enumerate some other discrete, microscopic models belonging to that class. We also discuss 
some curious properties of systems with complex noise, and analyze whether a real Langevin 
equation can be written for systems like the annihilation reaction A + A — > 0. 

I. ANALYSIS OF THE SINGLE-SITE EFFECTIVE POTENTIAL 

A. Reggeon field theory 

We start by analizing the zero-dimensional (single variable) version of the RFT Langevin 
equation. The Fokker-Planck equation associated to eq. ([l|) is [|14]: 



dP(n,t) d . , 2 ,„, . D9 2 n/ , 



By imposing the detailed balance condition, the associated formal stationary probability 
distribution is found to be : 

1 2 

P(n) = exp(— V(n)) oc — exp[— (an — bn 2 /2)] (5) 

where V(n) is the effective potential. In figure [I] we plot V(n) for different values of a, 
and D = 1; for o > y/2b the potential has a minimum at n ^ 0, while for a < -\/2& V(n) 
has no maximum or minimum. Note, however, that due to the (non-integrable) singularity 
at n = 0, the probability eq. (f|) is not normalizable, and the only stationary solution is 
P(n) = S(n). Therefore there is no active phase in this simple 0-dimensional case, and the 
systems decays towards the absorbing state for any set of paramenter values. 

Let us now study how the single-site effective potential behaves in dimensions larger than 
zero. In particular, we perform a numerical simulation of eq. ([I]) in one dimension. To do 



so we employ a technique developed by Dickman |25| to deal with numerical simulations of 
the continuous RFT. Let us point out that the simulation of this continuous theory with an 
absorbing state is not a trivial issue, and that, for example, a straightforward discretization 
of eq. (|5|) in which eventual negative values of the field (that may appear due to the 
discretization) are fixed to n(x) = 0, does not preserve the presence of an absorbing state. 
Dickman's method consists of a discretization of the space, time and also of the field variable, 



that ensures the presence of an absorbing state (see [25] for details). 

In the one-dimensional case the active phase survives to the effect of fluctuations contrar- 
ily to what happens in the single-variable case. In figure || we show the effective potential 
(defined as minus the logarithm of the normalized stationary probability distribution) for 
different values of a. The uppermost curve corresponds to a value of a in the absorbing 
phase; the second one to a = a critica i while the two lower ones are in the active phase. Note 
that in all the cases a singularity at the origin of the same type is present, and consequently, 
for any finite system there is a finite probability for the system in the active phase to go 
through the potential barrier and decay towards the absorbing state: the active phase is 
a metastable state. The mean time required for the system to overpass the barrier and 
collapse to the absorbing state grows exponentially with time, and becomes infinite in the 



thermodynamic limit. In this way the phase transition appears only in infinitely large sys- 
tems, and the large system-size limit has to be taken first than the infinite time limit in 
order to permit the presence of an active phase. In dimensions larger than d — 1 the same 
qualitative type of behavior is expected. 

B. Multiplicative noise 

The Fokker-Planck equation associated to eq. (Bt) in the zero- dimensional case is: 



,2 



dP(n.t) d , , 9 . „. . D d 

_Lli = __(„„ _ frn2)P(n , t) + -^n'P(n, t). (6) 

By imposing the detailed balance condition, the stationary formal solution is: 

1 — 2bn 

P{n) = exp(-y(n)) oc ^ 2(D ^ a)/D exp[-^-] (7) 

where V(n) is the effective potential; the solution is not normalizable when a < D/2, and 
normalizable otherwise. V(n) is plotted in figure ^|for different values of a and D — 1; the 
two lowermost curves correspond to eq. (0) in the absorbing phase (a = and a = 0.5) 
(where the only stationary solution is P(n) = 5(n)). The central one (a = 1), and the two 
uppermost curves (a = 1.5 and a = 2) are in the active phase. Note that contrarily to the 
RFT the MN exhibits a phase transition even in zero dimensions. Observe also that the 
singularity at the origin in the formal solution eq. (^) changes its degree as a is increased, 
in contrast with what happens for eq. (|5|); in fact, for D/2 < a < D (in the active phase) 
the singularity is integrable, and above a = D (also in the active phase) the origin becomes 
repelling instead of absorbing. This is an essential difference with RFT. 

Let us now explore how this property of the single-site potential is modified in higher 
dimensions. For that, we perform a numerical integration of the stochastic equation defining 



the model, which presents less technical difficulties than the integration of the RFT [22|. The 
result are presented in figure |j. Qualitatively the potential shape changes in the same way 
as it does in the 0-dimensional case. Above the critical point, there is either an integrable 
singularity (uppermost curve) or a repelling wall (three other curves) at the origin. The first 
case, i.e., an integrable singularity at the origin occurs in a very tiny region of the parameter 

6 



space. On the other hand, in the absorbing phase there is a collapse of the probability 
towards P(n) = 5(n) (non-integrable singularity at the origin). 

Therefore, the physics is very different from that of RFT: when the system is in the active 
phase, there is either an integrable singularity at the origin of the potential or a repelling 
wall. In any case, the situation differs from that in RFT, in which there are two locally 
stable attractors in the active phase: one at the origin (the stable one) and one at a different 
point (a metastable one). 

That is the reason why a reaction-diffusion system can not be described by an equation 
as eq. (|3|): in a finite reaction-diffusion system with a non- vanishing particle annihilation 
rate there is a non- zero probability of reaching the absorbing (empty) state for any finite 
system-size and for any set of parameter values. In systems with MN there is no accessible 
absorbing state in the active phase, i.e. there is no non-integrable singularity at the origin 
of the potential (and consequently no collapse of the probability density to the origin), and 
therefore MN does not capture the physics of reaction-diffusion systems. 

That fact, does not mean that it is not possible to construct discrete lattice model in 
the multiplicative noise universality class. In particular, systems exhibiting an unbinding 
transition from a wall (as, for example, the problem of local alignment of DNA chains 



and wetting transitions p6| ) belong in this universality class. These models are usually 
defined in terms of a field variable h(x,t) = ±log(n) that flows to =poo(n = 0) in the 
absorbing phase, and that reach a non-vanishing stationary average value otherwise. 

C. Peliti's field theory 

In this section we compare the effect of the complex noise appearing in the Peliti's field 
theory eq. (|2]), with the two previously studied cases to get a global picture of the different 
type of noises that can appear in systems with absorbing states. In order to understand 
what is the origin of the complex noise in processes as A + A — > 0, and to clarify whether 
a microscopic process like that with a real density field has necessarily to be described by 
a complex Langevin equation, we present a derivation of the Peliti's field theory for the 



annihilation reaction: A + A — * by employing the exact Poisson representation introduced 
by Gardiner and Chatuverdi p?| , pO"f . Note that for this reaction there is no active phase, 
and the interesting magnitudes are those describing the decay towards the absorbing state. 
For the sake of simplicity in the notation we present here the O-dimensional case, extensions 
to higher dimensions being straightforward. 
The master equation defining the process is: 

dP ^ t) = k[(n + 2)(n + l)P{n + 2, t) - n{n - l)P{n, t)} (8) 

Multiplying both sides of eq. (§) by s n , summing over all n's from to oo, and defining 
the generating function G(s,t) = J2™=o s n P(n, t), we get 

k(l-s 2 )d 2 s G(s,t). (9) 



9G(s,t) , ., 



dt 

We now introduce the Poisson transformation: 



F{ „. / , = J da QneXP [ Q h (a, t) (10) 



where /(a,t) is a given function (see 27[), in terms of which : 



G(s,t) = daf(a,t)exp(a(s — l)). (11) 

The Poisson transformation has the interesting property that the moments of P{n) and 
f(a) can be easily related: < a p >=< n(n— l)...(n— p+1) >', in particular the first moments 
are the same for both distributions. The integral over a can be taken over different domains 
of integration; for the moment let us assume a to be a real variable and leave the integration 
domain undetermined. In terms of f(a,t), eq. (|) reads, 

f daexp(a(s - l))d t f(a) = -k f daa 2 f(a,t)[-d 2 a + 2d a ] exp(a(s - 1)) (12) 

which, integrating by parts, and assuming that the boundary terms give a vanishing contri- 
bution to the integral |53| can be written as: 

d t f(a) = k[2d a a 2 f(a } t) - d 2 a a 2 f(a, t)] (13) 



which is a Fokker-Planck equation with a negative diffusion coefficient. The Langevin equa- 
tion stochastically equivalent to the previous Fokker-Planck equation is: 

d t a{t) = 2a{tf + iV2arj{t) (14) 

where rj(t) a Gaussian noise with amplitude 1, i the complex unit, and k has been eliminated 
by redefining the time as t — > t/k. Note that we have arrived to an inconsistency: a was 
assumed to be a real variable and we have arrived to a complex equation (observe that due 
to the complex term in eq. (|i"4|) a develops an imaginary part even if it is taken to be real 
at time t = 0. In appendix A we present further details on the impossibility of defining a 
real Poissonian representation for the reaction A + A — > 0). Let us now repeat the previous 
program but performing a complex transformation instead of a real one, i.e., we take J da 
to be /f^ da x Jf^ da y , where a x and a y are the real and imaginary parts of a respectively. 
This type of transformation leads to a function f(a) which is positive and can be identified 



as a probability distribution pTfl. Proceeding in that way we get a new set of Langevin 



equations for the variables a x and a 



u 



d t a x {t) = -2{a 2 x {t) - a 2 y {t)) + V2a y (t)r](t) 

d t a y (t) = -Aa x (t)a y (t) - V2a x (t) V (t) (15) 

that is equivalent to the original master equation. Note that both of the equations in eq. 
(|j~5l) include the same noise function r](t), and that they could be obtained straightforwardly 
from eq. (|i"4|) just by writing a = a x + ia y and separating the real and imaginary parts). 
A typical trajectory of the previous set of equations in the stationary state is shown in 
figure H it wonders in the complex plane avoiding a region around 0. Even if the stationary 
solution of the underlying process A + A — > is a delta function at zero or one particles 
(depending on whether the initial condition is even or odd respectively), the stationary 



probability associated to eq.fll^) is not a delta function, but some complicate distribution 
with < a x >= 1/2 and < a y >= 0. The value 1/2 comes from the fact that for initial 
conditions with n even, n(t — > oo) = and, for n odd, n(t — > oo) = 1; the variable 
< a x >=< n > is the average of the two previous possibilities. On the other hand, the 
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expectation value of the imaginary part is zero as expected given the relation among moments 
of f{a) and P{n). It is interesting to note that the effective potential associated to the 
stationary distribution is a non-differentiable one; in figure (|^) we show a one-dimensional 
cut of the stationary potential for different a x with a y = as computed in a simulation of 
eq. flI5p. Note that contrarily to the cases of the RFT and multiplicative noise equation, 
now the dynamics is not frozen even if the system has relaxed to the absorbing state. 

We now explore the possibility of finding a real-variable Langevin equation describing 
this class of systems |29| . 



As the second equation in eq. (|T^) is linear in a y it is possible to integrate it analytically; 
doing so and substituting the result in the first one, we get a closed equation for a x , that 
reads: 

d t oc x (t) = a x - 2a 2 x (t) + 21{tf + V2I(t)r](t) (16) 

with 

I(t) = a y (0) exp(- f dt\Aa x - 1)) - V2 f dt'a x {t')r]{t') exp(- f dt"(Aa x - 1)) (17) 

which is a non-Markovian equation (see appendix B). The stationary potential associated 
to eq.(|16|) cannot be calculated analytically; the numerical solution is shown in figure [7]. 
First we observe that it is non-differentiable at a x = 1/2; also we point out that a x is not 
absorbing in general (except for the pathological and unphysical case a J/ (0) = 0), in other 
words: due to the presence of the non-Markovian terms, proportional to I(t), the system 
can cross from positive values to negative ones. 

The role of the complex variable in eq. flUp is played, after a y has been integrated 
out, by the non-Markovian terms in eq. (0), and in both cases the absorbing state of 
the microscopic associated process is not described by a frozen dynamics in the Langevin 
representation, but by a non-trivial dynamics (complex or non-Markovian) which statistical 
properties reproduce those of the reaction-diffusion model. Therefore, the nature of the 
absorbing state in this case is essentially different from those of the previously studied cases. 

For the sake of completeness let us just mention briefly that a numerical study of the one- 
dimensional Peliti's field theory in terms of a complex Langevin equation has been recently 

10 



published ||30|| . The measured magnitudes are in very good agreement with the theoretical 



predictions coming from renormalizat ion-group and other type of analysis p!8| , pl| , confirming 
that a complex representation captures the physics of microscopic systems as, for example, 
the procces A + A — > 0. 

II. CONCLUSIONS 

We have analyzed different Langevin equations associated to systems with absorbing 
states. Systems described by the Reggeon field theory Langevin equation exhibit a non- 
integrable singularity at the origin of the single-site potential, that corresponds to a true 
absorbing state, i.e., there is an accumulation of probability density at the origin, while 
the active state is a metastable one for finite system-sizes. Systems with multiplicative 
noise instead change the degree of the singularity at the origin as the control parameter is 
changed: while in the absorbing phase there is a collapse of the probability density towards 
the origin, in the active phase there is either an integrable singularity at the origin or 
it becomes repelling, in which case the the probability to be nearby the origin becomes 
extremely small, and there is no accessible absorbing state. That is the reason why is not 
possible to find reaction-diffusion systems (in which for finite size systems there is always 
a finite probability of reaching the absorbing state) in the multiplicative noise universality 
class. We have also analized some aspects of the annihilation process A + A — > 0, which is 
described by a complex noise Langevin equation or alternatively by a real non-Markovian 
Langevin equation. This type of Langevin equation shows a behavior quite different from 
that of Reggeon field theory and multiplicative noise; in particular, even if the system is in 
the absorbing state, there is not a collapse of the probability density to a delta function, 
and the dynamics is non-trivial. Systems with complex noise can alternatively be described 
by real non-Markovian equations. 
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III. APPENDIX A 

Let us consider the pair of reactions A + A — > and A — > 2A, the first occurring with a 
rate k 2 and the second with k\. The associated master equation is 

dP(n,t) 



dt 



k 2 [(n + 2)(ra + l)P(n + 2, t) - n(n - l)P(n, £)] 



+ fci[(n- l)P(n-l) -nP(n)]. (18) 

Performing a real Poissonian transformation we get: 

dtf(a) = [d a (ha - 2k 2 a 2 )f(a, t) + d 2 a (k x a - k 2 a 2 )f(a, £)] (19) 

which is equivalent to the Langevin equation 

d t a{t) = (ha - 2k 2 a 2 ) + V2(k ia - k 2 a 2 ) 1/2 r](t) (20) 

interpreted in the Ito sense. Note that the factor multiplying the noise is positive in the in- 
terval a g]0, ki/k 2 [, and vanishes at the limits of the previous interval. The formal stationary 
solution of eq. (|19"D is : 

f(a) oc - exp(2a)(l - affi-W (21) 

a 

with j3 = k 2 jk\. Considering the Poisson representation as defined in [0, ki/k 2 ], it is a matter 
of simple algebra to verify that the boundary terms appearing in the processes of getting 
eq. (|T9|) from eq. fll8|) give a vanishing contribution. At the same time, trajectories of eq. 



(20|) with initial condition in [0, k\jk 2 \ do not leave that interval. On the other hand, if the 
domain of integration was extended over those limits, eq. (|2~Cf) would develop an imaginary 
part and the procedure would not be self-consistent. Therefore the transformation is well 
defined only in the real interval [0, k\jk 2 \. From a renormalization group point of view the 
noise-term proportional to k 2 in eq. ( pD| ) can be argued to be irrelevant rendering the system 
in the RFT universality class. 

We can now take the limit fci — > to see what happens in the Peliti's field theory case: 
the interval in which the Poisson representation is defined shrinks down to a single point; 
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a = 0. In the strict limit k\ = 0, a meaningful real Poisson representation can not be 
performed, and a complex representation is required. 

Let us point out as a final remark that it is somehow surprising that the standard renor- 
malization group analysis of the Peliti's field theory, based on a path integral representation 
of eq. (13) (or equivalently of eq.flTBD), in which a is treated as a real variable, give the right 



exponents and properties |L8],[[l|. We will investigate that apparent paradox in a future 
work. 

IV. APPENDIX B 

As a last attempt to write down a one- variable Langevin equation with a structure simpler 
than eq. (0), and inspired by the rotational quasi-symmetry of the stationary distribution 
solution (see figure |5|), we perform a change of variables to polar coordinates p and 6 defined 
by: a x = pcos(6), a y = psin(9). After changing variables (for which Ito calculus is required 
0), we get: 

d t p(t) = p(t) - 2p 2 (t) cos(9) 

d t 9{t) = -2p(t) sin(0) - V2 V {t). (22) 

Observe that the first one is a deterministic equation, while the second one is stochastic. It 
is easy to verify that this system does not admit a potential solution (which is consistent 
with the stationary potential being non-differentiable). This new set of equations permits to 
derive analytically some of the properties of the stationary probability distribution (as, for 
example, the presence of a maximum at (1/2,0)), but it does not simplify the elimination 
of one of the variables in favor of the other one to construct a simple one-variable Langevin 
equation. 
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FIGURES 
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FIG. 1. Potential V{n) associated to the formal stationary solution of the O-dimensional RFT 
for different parameter values: b = D = 1, and from top to bottom: a = 0.4,0.8,1,1.5 and 1.9. 
Note the presence of a strong (non-integrable) singularity at the origin in any case. 




4.0 



FIG. 2. Stationary potential for the one-dimensional RFT as coming from a simulation of the 
discretized Langevin equation. The uppermost curve corresponds to a value of a in the absorbing 
phase; the second one to a = a C ritical while the two lower ones are in the active phase. 
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FIG. 3. Potential V(n) associated to the formal stationary solution of the O-dimensional 
multiplicative noise equation for different parameter values: b = D = 1, and from top to bottom: 
a = 2,1.5,1,0.5 and 0. The potential develops a minimum as a is increased, has a negative 
singularity at the origin for a < 1, and a positive singularity for a > 1; the singularity is integrable 
in the active phase, i.e., when a < 0.5, while in the absorbing phase the only stationary solution is 
a delta function at the origin. 
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FIG. 4. Stationary potential for the one-dimensional multiplicative noise Langevin equation 
as coming from a simulation of the discretized equation. The curves correspond to four different 
parameter values all of them in the active phase. Note that the singularity at the origin is positive 
for the three lower most curves (a = —2, —2.1 and —2.2 respectively), therefore there is not an 
absorbing, but a repelling state. The uppermost curve (a = —2.23), with a negative singularity at 
the origin, is still in the active phase, but the singularity is integrable. 
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FIG. 5. Typical trajectory of the complex noise equation in the stationary regime. Note that 
there is a large probability of finding the system in the vicinity of 1/2, 0, while points inside a circle 
centered at the origin of radius 1/2 are inaccessible. 
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FIG. 6. Cut of the stationary probability distribution associated to the complex noise equation 
with a y = 0. Observe the non-differentiability at a x = 1/2. 
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FIG. 7. Stationary probability distribution function associated to the non-Mar kovian equation 
or equivalently, the projection of the stationary probability function associated to the complex noise 
equation over the real axis. 
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